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Abstract 

A  coupled  wave-hydrodynamic  modeling  system,  comprised  of  a  random 
wave  model  driving  a  quasi-3D  nearshore  hydrodynamic  model,  is  described. 
Random  wave  formulations  for  several  inputs  to  the  hydrodynamic  model  arc 
developed.  An  alternate  wave  dissipation  mechanism  is  incorporated  into  the 
random  wave  model,  and  two  wave  roller  descriptions  arc  implemented  to  cal¬ 
culate  volume  flux  and  other  roller-dependent  input  properties.  Comparison  to 
laboratory  and  field  data  indicate  that  an  evolving  roller  description,  in  conjunc¬ 
tion  with  the  3D  dispersive  mixing  inherent  in  the  hydrodynamic  model,  yield  the 
best  results.  A  method  to  nest  the  model  system  inside  larger-scale  wave  models 
is  described,  and  an  application  to  an  area  of  complex  bathymetry  shown. 


INTRODUCTION 

Numerical  modeling  of  wave-induced  nearshore  circulation  has  undergone  much 
development  in  recent  years,  to  the  extent  that  use  in  field  situations  is  almost  rou¬ 
tine.  The  one-dimensional  Navy  Standard  Surf  Model  (Earle  1989)  has  been  used  for 
forecasting  nearshore  conditions  for  Navy  operations  for  several  years.  The  Delft3D 
model  (http://www.wldelft.nl/soft/d3d)  is  a  commercial  package  used  in  engineering 
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applications  on  a  worldwide  basis.  More  recently,  the  quasi-3D  model  SHORECIRC 
(Van  Dongeren  and  Svendsen  2000)  has  been  validated  with  data  from  the  Duck94  ex¬ 
periment  (Svendsen  et  al.  1997).  It  is  therefore  opportune  to  investigate  incorporation 
of  wave  forcing  more  in  line  with  general  nearshore  conditions  in  the  field. 

In  this  paper  we  discuss  the  incorporation  of  random  wave  effects  into  the  quasi-3D 
model  SHORECIRC.  Required  wave  forcing  for  the  model  extends  past  the  standard 
gradients  of  radiation  stress;  random  wave  extensions  of  these  inputs  are  developed 
herein.  Comparisons  to  data  (laboratory  and  field)  are  shown,  and  the  incorporation  of 
the  revised  model  into  a  model  nesting  scheme  is  described. 

SHORECIRC  and  REF/DIF-S 

The  SHORECIRC  model  (Van  Dongeren  and  Svendsen  2000)  combines  the  phys¬ 
ical  effects  of  a  non- depth- uniform  current  structure  with  a  two-dimensional  model¬ 
ing  formalism,  thereby  requiring  only  horizontal  discretization  of  the  domain.  Semi- 
analytic  solutions  are  used  to  replicate  the  depth  structure  of  the  currents,  which  is  then 
integrated  to  provide  the  necessary  coefficients  for  the  current-current  and  current- 
wave  interaction  terms  in  the  two-dimensional  horizontal  model.  This  allows  for  the 
effect  of  3D  dispersive  mixing  (Putrevu  and  Svendsen  1999),  essential  for  obtaining  the 
proper  horizonal  momentum  exchange  in  the  nearshore  without  requiring  unrealistic 
values  of  eddy  viscosity  coefficients  for  sufficient  lateral  mixing. 

The  scales  of  motion  simulated  by  SHORECIRC  are  assumed  to  be  averaged  over 
individual  waves  (wave  group  scale)  and  thus  wave  information  is  required  from  a 
separate  model.  This  information  includes  gradients  of  radiation  stress;  short  wave 
volume  flux;  near-bed  orbital  velocities;  and  estimates  of  breaking  wave  dissipation. 
Short  wave  volume  flux  is  required  due  to  the  quasi-3D  nature  of  the  model  and  the 
need  to  balance  the  shoreward  mass  flux  with  the  resulting  undertow. 

The  default  version  of  the  SHORECIRC  model  presently  available  uses  the  single 
frequency  short  wave  model  REF/DIF- 1  (Kirby  and  Dalrymple  1994)  as  a  wave  driver. 
The  REF/DIF- 1  model  uses  the  parabolic  approximation  to  the  mild-slope  equation 
(Berkhoff  1972;  Radder  1979)  for  wave  propagation  over  varying  bathymetry,  and  is 
thus  capable  of  simulation  both  refraction  and  diffraction,  with  higher  order  terms  in¬ 
cluded  to  enhance  model  accuracy  at  oblique  angles  (Kirby  1986).  While  useful  for 
development  purposes  and  research  applications,  general  utility  is  hampered  with  a 
monochromatic  wave  driver.  The  model  REF/DIF-S  (Chawla  et  al.  1998)  is  a  spectral 
version  of  REF/DIF- 1;  it  propagates  many  frequency  and  direction  components  over 
the  domain  simultaneously  within  the  REF/DIF- 1  modeling  format,  using  the  dissipa¬ 
tion  function  of  Thornton  and  Guza  (1983)  to  represent  energy  decay  in  the  spectrum 
due  to  breaking.  Both  REF/DIF- 1  and  REF/DIF-S  are  wave-resolving  models,  and 
thus  directional  spectra  must  be  decomposed  into  individual  wavetrains  with  random 
phases  for  initialization.  For  later  reference,  we  write  the  REF/DIF-S  governing  equa¬ 
tion  in  generic  form  as: 

CgnAn  +  MAnyy  +  N An  -f  P Anxy  +  QAnyyx  +  aAn  =  0  (1) 
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where  Cgn  is  the  group  velocity  for  wave  component  n;  M,  N,  P  and  Q  are  coef¬ 
ficients  for  various  transformation  processes  representing  wave  shoaling,  refraction, 
diffraction,  wave-current  interaction  and  wide-angle  propagation  effects;  An  is  the 
complex  amplitude  of  the  n’th  wave  component  in  the  spectrum;  and  a  is  the  dissi¬ 
pation  function  used  to  simulate  wave  energy  decay  due  to  breaking.  Subscripts  x,  y 
refer  to  partial  differentiation. 

WAVE  FORCING 

The  REF/DIF- 1  model  used  in  the  default  version  of  SHORECIRC  provides  the 
hydrodynamic  model  with  all  necessary  forcing  components.  It  uses  the  results  of 
Svendsen  (1984)  to  calculate  short  wave  volume  flux  inside  and  outside  the  surf  zone, 
and  uses  Dally  et  al.  (1985)  to  provide  dissipation  estimates.  Since  REF/DIF-1  is 
monochromatic,  a  demarcation  between  non-breaking  and  breaking  regions  is  possi¬ 
ble,  and  surf  zone  mechanisms  can  be  switched  on  or  off  accordingly.  This  division 
is  not  possible  in  irregular  wave  models  such  as  REF/DIF-S.  While  the  probabilistic 
dissipation  function  of  Thornton  and  Guza  (1983),  included  in  the  default  version  of 
REF/DIF-S,  works  well  for  simulating  waveheight  decay,  there  is  no  clear  methodol¬ 
ogy  for  incorporating  or  excluding  surf  zone  processes  due  to  this  lack  of  a  “break¬ 
point;”  this  problem  needs  to  be  addressed.  Additionally,  random  wave  equivalents  for 
several  input  quantities  such  as  short  wave  volume  flux  require  development. 


Dissipation  Mechanisms 

The  default  dissipation  mechanism  in  REF/DIF-S  is  that  of  Thornton  and  Guza 
(1983).  This  function  was  adapted  for  complex  amplitude  models  by  Chawla  et  al. 
(1998): 


=  30F  fpB 3  5 
4  7  4/i5  rms 


(2) 


where  B  and  7  are  free  parameters  and  fp  is  the  peak  frequency.  The  dissipation 
function  of  Battjes  and  Janssen  (1978)  was  added  to  REF/DIF-S  as  an  alternative.  It 
was  adapted  for  complex  amplitude  models  by  Eldeberky  and  Battjes  (1996): 


1  fpQbHmaX 

4£CsnjAn|2 


(3) 


where  Qb  is  the  percentage  of  breaking  waves,  calculated  by  the  following  implicit 
relationship: 


/ 1  —  Qb\  _  _  f  Hrms  \  2 
\  In  Qb  )  V  Hmax  ) 


Hr, 


0.88 

T 


tanh 


(5) 


where  7'  is  a  parameter  dependent  on  offshore  wave  steepness,  and  k  is  the  wavenum¬ 
ber  associated  with  the  characteristic  frequency  for  the  spectrum  (usually  the  peak 
frequency).  Either  dissipation  function  is  a  selectable  option  in  the  model. 
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Roller  Descriptions 

In  order  to  calculate  the  forcing  and  short  wave  volume  fluxes,  we  will  require  ex¬ 
pressions  for  the  roller  area  (per  unit  length  of  crest)  A'  and  the  roller  energy  density 
Er.  These  are  in  turn  dependent  on  the  particular  description  of  the  roller  characteris¬ 
tics  selected. 


In  REF/DIF-S  we  make  use  of  two  roller  descriptions.  The  first  is  denoted  the 
“static  roller”  because  it  is  dependent  entirely  on  local  water  depth  properties.  We  use 
the  result  of  Fippmann  et  al.  (1996),  who  extended  the  roller  area  expression  for  weak 
hydraulic  jumps  to  include  random  wave  effects  by  integrating  the  cube  of  the  breaking 
waveheight  through  the  probability  density  function  for  breaking  waves  (Thornton  and 
Guza  1983).  The  resultant  expression  for  roller  area  A'  is: 


where: 


A,  ( 

4/rtana 


/ 

/fj3\  _  3^  1  rj5 

'  b/  4  (jh)2  rms 

V 


\ 

/ 


(6) 


(7) 


and  a  is  the  angle  of  the  stress  vector  at  the  wave/roller  interface.  The  roller  energy 
density  Er  is  found  from  (Svendsen  1984): 


Er 


pA'C 
2  L 


(8) 


where  C  is  the  average  phase  speed  and  L  is  the  average  wavelength. 


The  second  description  uses  the  evolving  roller  model  of  Stive  and  deVriend  (1994). 
The  governing  equation  for  the  evolving  roller  is: 

d(2ErCcos9)  —  2gEr  sin 
dx  C 

where  S  is  the  dissipation  from  wave  breaking  (using  the  original  form  of  dissipation 
prior  to  adaptation  for  REF/DIF-S)  and  fi  is  the  slope  of  the  roller  interface.  This  model 
is  solved  using  a  centered  finite  difference  scheme  solved  at  every  step  of  the  REF/DIF- 
S  model  to  yield  Er  through  the  domain.  The  model  is  initialized  with  Er  =  0  at  the 
offshore  boundary.  The  roller  area  A!  is  then  calculated  from  Er  using  (8). 

Radiation  Stress  in  the  Surf  Zone 

The  radiation  stress  evaluations  in  REF/DIF-S  essentially  sum  contributions  to  the 
overall  radiation  stress  from  each  frequency  and  direction  component  in  the  spectrum. 
Inside  the  surf  zone,  however,  a  randomization  of  the  expression  of  Svendsen  (1984) 
for  radiation  stress  due  to  rollers  is  used: 


Sxxr  =  2 Er  cos2  9 


(10) 
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Syyr  =  2  Er  sin2  9 


(ID 


Sxyr  =  2 Er  cos  9  sin  9  (12) 

where  9  is  the  mean  angle.  This  is  then  added  to  the  radiation  stress  contributions  from 
the  non-breaking  waves. 

Short  Wave  Volume  Flux 

Outside  the  surf  zone,  short  wave  volume  flux  is  essentially  the  Stokes  drift.  Within 
the  surf  zone,  the  form  of  the  volume  flux  depends  on  the  assumed  characteristics  of 
the  breaking  wave.  For  monochromatic  waves,  Svendsen  (1984)  determined  a  volume 
flux  assuming  that  the  breaking  wave  can  be  described  by  a  “roller”  of  white  water 
rotating  on  the  face  of  an  unbroken  wave.  For  random  waves,  we  assume  that  the  short 
wave  volume  flux  can  be  described  generically  as: 

A'  - 

Qwrx  —  -=  COS  9  (13) 

A!  - 

Qwry=YSin0  (14) 

where  T  is  the  average  wave  period.  This  is  then  added  to  the  Stokes  drift  to  obtain  the 
total  short  wave  volume  flux. 

Eddy  viscosity  specification 

As  mentioned  previously,  the  quasi-3D  nature  of  the  SHORECIRC  allows  proper 
mixing  to  take  place  without  requiring  egregiously  large  values  of  eddy  viscosity  coef¬ 
ficient.  Thus  the  turbulent  shear  stresses  contribute  litte  to  the  overall  mixing.  The  eddy 
viscosity  model  used  for  the  turbulent  shear  stresses  is  depth  uniform;  the  coefficient 
is: 


vt  =  C1K^^u0h  + Mh  +vto  +  vs  (15) 

where  k  is  von  Karman’s  constant  0.4);  C\  and  M  are  coefficients  for  bottom- 
induced  turbulence  and  breaking  turbulence,  respectively;  fw  is  the  wave-related  bot¬ 
tom  friction  coefficient;  ua  is  the  magnitude  of  the  near-bed  particle  velocity;  and  vto 
and  vs  are  the  background  and  Smagorinsky  subgrid  viscosities,  respectively.  While 
the  first  term  of  (15)  is  relevant  for  the  entire  domain,  the  second  term  only  pertains 
to  the  surf  zone.  Since  the  dissipation  mechanism  is  always  activated  in  REF/DIF-S, 
all  terms  in  (15)  are  active  for  the  entire  domain,  with  the  dissipation  a  weighting  the 
contribution  of  the  breaking  wave  turbulence.  Additionally  an  estimate  of  ui/3  replaces 
u0  in  REF/DIF-S. 

COMPARISON  TO  DATA 

In  this  section  we  compare  the  results  of  SHORECIRC  from  the  different  forcing 
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mechanisms  to  data.  Our  first  comparison  is  to  the  data  of  Reniers  and  Battjes  (1997), 
who  ran  both  monochromatic  and  irregular  waves  over  a  sloping  bathymetry  with  a 
nearshore  bar.  Waveheights,  mean  sea  surface  and  nearshore  depth-averaged  currents 
were  measured.  They  compared  their  results  to  a  one-dimensional  wave-current  model, 
with  good  results. 

Here  we  compare  the  SHORECIRC  model  driven  with  both  the  Thornton  and  Guza 
(1983)  and  Battjes  and  Janssen  (1978)  dissipation  models.  Both  the  static  roller  and 
the  Stive  and  deVriend  (1994)  roller  are  used.  Figure  1  shows  a  comparison  of  mod¬ 
eled  currents  to  data.  The  modeled  currents  resulted  from  the  following  combinations: 
Battjes  and  Janssen  dissipation  with  Stive  and  deVriend  roller;  Thornton  and  Guza  dis¬ 
sipation  with  Stive  and  deVriend  roller;  and  Thronton  and  Guza  dissipation  with  static 
roller.  All  tests  were  performed  with  the  quasi-3D  dispersive  mixing  mechanism  acti¬ 
vated.  Free  parameters  for  the  Thornton  and  Guza  decay  were  B  =  0.8  and  7  =  0.6. 
For  the  static  roller  a  =  10°,  while  for  the  evolving  roller  0  =  0.1.  The  bottom  friction 
coefficient  fw  =  0.02,  while  the  coefficients  for  (15)  are  C^k  =  0.08  and  M  =  0.08. 
Both  dissipation  mechanisms  performed  well  in  capturing  the  waveheight  decay,  with 
the  Thornton  and  Guza  result  exhibiting  a  slightly  better  comparison.  However,  the 
Battjes  and  Janssen  decay  with  the  Stive  and  deVriend  roller  yielded  the  most  favor¬ 
able  comparison  to  velocity  data;  the  Thornton  and  Guza  decay  with  the  Stive  and 
deVriend  roller  also  compared  reasonably  well.  On  the  other  hand,  the  static  roller 
resulted  in  an  unrealistically  narrow  longshore  current  profile;  it  is  apparent  that  the 
evolving  roller  of  Stive  and  deVriend  (1994)  yields  a  more  realistic  longshore  current 
profile.  Tests  with  the  quasi-3D  mechanism  deactivated  (not  shown)  clearly  exhibited 
insufficient  mixing,  even  with  increased  values  of  eddy  viscosity  coefficient. 

Comparisons  to  field  data  from  the  Duck94  experiment  (held  at  Duck,  NC  during 
October  1994)  were  performed  to  evaluate  the  model  under  field  conditions.  Nearshore 
waveheights,  undertow  profiles  and  depth  averaged  currents  were  compared  to  mea¬ 
surements  for  the  case  of  October  12,  1994.  Both  a  high-tide  condition  (12:30  EST) 
and  a  mid-tide  condition  (16:20  EST)  were  simulated  with  the  REF/DIF-S-forced 
SHORECIRC  model,  using  the  Battjes  and  Janssen  decay  with  the  Stive  and  deVriend 
roller.  Figures  2  and  3  show  the  high  tide  case.  From  Figure  2  it  is  apparent  that  the 
depth  profiles  for  the  longshore  currents  are  well  simulated.  The  cross-shore  current 
comparison  does  evidence  some  deviation  from  measurements  over  the  bar  near  the 
surface.  Figure  3  shows  comparisons  to  the  nearshore  waveheights  and  depth- averaged 
longshore  currents.  The  model  appears  to  underpredict  the  amount  of  dissipation  in  the 
wavefield.  This  insufficient  dissipation  over  the  bar  is  coincident  with  the  underpredic¬ 
tion  of  the  undertow  curvature  and  the  maximum  longshore  current.  This  is  presently 
being  investigated;  however,  the  level  of  agreement  displayed  here  is  encouraging. 

MODEL  NESTING 

To  fiurther  increase  the  general  utility  of  the  REF-DIF/S-SHORECIRC  modeling 
system,  a  provision  for  accepting  initial  and  boundary  conditions  from  offshore  models 
is  included.  Generally,  initial  conditions  for  nearshore  wave  models  are  provided  by 
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0.08 


Figure  1:  Comparison  to  experiment  of  Reniers  and  Battjes  (1997).  Top:  Hrms ;  solid 
line:  Battjes  and  Janssen  decay;  dashed  line:  Thornton  and  Guza  decay.  Bottom: 
Longshore  current;  solid  line:  Battjes  and  Janssen  decay  with  Stive  and  deVriend 
roller;  dashed  line:  Thornton  and  Guza  decay  with  Stive  and  deVriend  roller;  dashed- 
dot  line:  Thornton  and  Guza  decay  with  static  roller.  Symbols  indicate  measurements. 
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Cross-shore  current 


Long-shore  current 


Figure  2:  Comparison  to  Duck94  undertow  measurements  for  October  12,  1994, 
1230EST  (high  tide).  Top:  Cross-shore  velocities.  Bottom:  Longshore  velocities. 
Solid  line:  results  from  REF/DIF-S-forced  SHORECIRC  using  Battjes  and  Janssen 
wave  decay  with  Stive  and  deVriend  roller.  Crosses:  data. 
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Battjes  and  Janssen/Stive  and  deVriend  -  DUCK94  data  (hightide:1 2:30,  Oct  12,  1994) 


Figure  3:  Comparison  to  Duck94  waveheight  and  longshore  current  measurements  for 
October  12,  1994,  1230EST  (high  tide).  Top:  Hrms.  Middle:  Depth  averaged  velocity. 
Bottom:  Depth  profile.  Solid  line:  results  from  REF/DIF-S-forced  SHORECIRC  using 
Battjes  and  Janssen  decay  with  Stive  and  deVriend  roller.  Asterisks:  waveheight  data. 
Open  circles:  depth- averaged  current  data. 
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large  scale  wave  propagation  models  such  as  SWAN  (Booij  et  al.  1999).  These  models 
are  typically  phase-averaged,  and  thus  some  transformation  of  variables  is  required 
for  input  to  the  phase-resolving  wave  model  REF/DIF-S.  This  problem  becomes  more 
involved  if  wave  conditions  along  the  offshore  boundary  of  the  nearshore  grid  change 
substantially. 

Since  the  REF/DIF-S  model  is  linear,  the  phase  relationships  between  different 
spectral  components  are  arbitrary;  however,  the  spatial  gradients  of  the  phase  for  each 
component  are  essential  since  they  relate  to  the  wave  angle.  Given  directional  spec¬ 
tra  at  every  point  along  the  offshore  boundary  of  the  nearshore  grid  S(f,  9)(y),  the 
transformation  to  complex  ampitudes  An  is: 

An{y)  =  ^2S(f,9)(y)AfAyyJx^dy  (16) 

where: 

Hf,6,y)  =  k(f,y)sm6  (17) 

and  where  the  maximum  value  of  n  is  the  number  of  frequencies  times  the  number  of 
directions.  This  has  the  effect  of  phase-lagging  the  offshore  wave  condition  for  each 
frequency  and  direction  component  along  the  offshore  boundary  of  the  nearshore  grid. 
In  this  manner  the  variation  due  to  wave  propagation  offshore  of  the  nearshore  grid  is 
properly  reflected  in  the  initial  condition. 

The  modeling  system  is  applied  to  the  bathymetry  of  the  upcoming  Nearshore 
Canyon  Experiment  (NCEX),  slated  to  begin  in  fall  2003  near  Fa  Jolla,  California.  The 
domain  and  waveheights  from  the  SWAN  model  are  shown  in  Figure  4;  it  is  apparent 
that  there  is  considerable  variation  in  the  nearshore  wave  climate.  The  bathymetry  of 
the  domain  was  gleaned  from  the  National  Ocean  Survey  (NOS)  at  3  second  resolu¬ 
tion.  The  resolution  of  the  SWAN  simulation  was  Ax  =  38.5m  and  Ay  =  46.5 m,  with 
the  nearshore  grid  (shown  in  the  small  box  in  Figure  4)  set  at  16  times  this  resolution 
(Ax  =  2.40625m,  Ay  =  2.90625m).  Figure  5  shows  the  nearshore  waveheights  and 
depth- averaged  circulation  from  the  REF/DIF-S  -  SHORECIRC  system.  Though  the 
bathymetry  in  this  area  is  quite  planar,  severe  bathymetric  variations  offshore  of  this 
area  exacts  a  large  degree  of  variation  on  the  waveheight  field,  resulting  in  a  circulation 
field  more  complex  than  the  planar  bathymetry  would  imply. 

CONCLUSIONS 

A  coupled  wave-hydrodynamic  model  system  was  described.  This  system  is  com¬ 
prised  of  a  random  nearshore  wave  model  (REF/DIF-S)  coupled  with  a  quasi-3D 
nearshore  hydrodynamic  model  (SHORECIRC).  Certain  enhancements  were  made  to 
the  REF/DIF-S  model  to  make  it  a  more  general,  flexible  model.  To  compliment  the 
existing  default  wave  dissipation  mechanism  of  Thornton  and  Guza  (1983),  the  wave- 
height  decay  formulation  of  Battjes  and  Janssen  (1978)  was  included  as  an  option. 
Random-wave  formulations  of  roller-induced  short  wave  volume  flux  and  radiation 
stresses  were  developed  and  implemented  into  the  REF/DIF-S  model.  These  quanti¬ 
ties  are  dependent  on  the  wave  roller  description  used.  We  detailed  two  mechanisms: 
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Figure  4:  Bathymetry  and  Hi/ 3  over  NCEX  bathymetry,  Scripps  Canyon,  CA. 
Nearshore  grid  outlined  on  right  side  of  figure  near  coast. 


Velocity  vectors 


Figure  5:  Waves  and  currents  over  nearshore  grid,  Scripps  Canyon.  Left:  H i/3.  Right: 
Current  vectors;  the  longest  vector  corresponds  to  1  m/s. 
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a  “static”  roller  whose  properties  are  only  based  on  local  conditions,  and  the  evolving 
roller  of  Stive  and  deVriend  (1994).  Additionally,  wave  energy  dissipation  was  used 
to  weight  the  surf-zone  contribution  to  the  overall  eddy  viscosity  coefficient  for  use  in 
modeling  turbulent  shear  stresses. 

Comparisons  to  laboratory  data  (Reniers  and  Battjes  1997)  and  field  data  from 
the  Duck94  experiment  revealed  that  evolving  roller  of  Stive  and  deVriend  (1994), 
along  with  the  quasi-3D  dispersive  mixing  in  SHORECIRC,  yields  the  best  results. 
Deactivating  the  dispersive  mixing  leads  to  insufficient  mixing  and  an  unrealistically- 
narrow  longshore  current  profile;  use  of  the  static  roller  yields  the  same  results. 

Finally,  a  capability  for  nesting  the  coupled  modeling  system  inside  a  larger  scale 
model  was  described.  A  method  for  transforming  directional  spectra  saved  along  the 
offshore  boundary  of  a  nearshore  grid  to  complex  amplitudes  was  detailed.  The  sys¬ 
tem  was  applied  to  the  site  of  the  upcoming  NCEX  experiment;  the  offshore  wave 
conditions  were  provided  over  the  large  scale  region  by  the  SWAN  model  (Booij  et  al. 
1999).  Directional  spectra  were  saved  along  the  offshore  boundary  of  the  nearshore 
grid  and  transformed  into  complex  amplitudes  for  the  REF/DIF-S  model,  the  results 
of  which  were  subsequently  input  to  SHORECIRC.  The  resulting  waveheights  and 
currents  over  the  nearshore  domain  showed  considerable  variation,  despite  the  rather 
planar  bathymetry  in  the  nearshore. 
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